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Abstract. Dynamical relativistic effects are often neglected in the description of reactions with unstable nuclear beams at 
intermediate energies (Ei a i, ~ 100 MeV/nucleon). Evidently, this introduces sizable errors in experimental analysis and theo- 
retical descriptions of these reactions. This is particularly important for the experiments held in GANIL/France, MSU/USA, 
RIKEN/Japan and GSI/Germany. I review a few examples where relativistic effects have been studied in nucleus-nucleus 
scattering at intermediate energies. 
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1. INTRODUCTION 

The number of radioactive beam facilities are growing fast around the world. Some of these facilities use the 
fragmentation technique, with secondary beams in the energy range Eutb ~ 100 MeV/nucleon. Examples are the 
facilities in GANIL/France, MSU/USA, RIKEN/Japan and GSI/Germany. Relativity, an obviously important physics 
concept [1], is often neglected in calculations aiming at relating the reaction mechanisms to the internal structure of 
the projectiles. For example, popular DWBA codes (FRESCO, ECIS, DWUCK, etc.) useful in the analysis of nuclear 
reactions, include relativity only in kinematic relations. The effects of relativity in the reaction dynamics (i.e. the 
interaction) is not accounted for because these codes were intended for lower energies. It is also important to notice that 
the inclusion of relativistic effects in the nucleus-nucleus dynamics is a very difficult task. A fully covariant treatment 
of the nuclear many-body scattering (with inclusion of retardation effects between all nucleons) is not possible without 
approximations. 

In this short article I will review a few examples where relativistic effects have been included in nuclear reactions 
at intermediate energies. It is worthwhile to observe that after 100 years of relativity [1] this still remains a challenge 
in many aspects. 



2. SEMICLASSICAL METHODS AND ELASTIC SCATTERING 

Semiclassical methods are a very popular tool for the description of nucleus-nucleus collisions at high energies. As an 
example, I cite the Coulomb excitation mechanism, in which the inelastic cross section can be factorized as 

where P,^f is the probability for a nuclear transition between the states i and / when the nuclei scatter through an angle 
Q.. The collision dynamics enters here in two distinct ways: in the calculation of Pi-^f and (da /dQ.) el . Let us forget 
about Pi-vf for a moment and let us investigate (da /dQ) A . The question here is what is the error done by measuring 
the number of particles scattered to Q. and using the Rutherford formula for (da/dD.) e] . This has been investigated in 
ref. [2]. 

A system of two point charges interacting electromagnetically and moving at low velocities can be described by 
an approximate Lagrangian which depends only on the degrees of freedom of the particles neglecting those related 
to the electromagnetic field (the Darwin Lagrangian's). In this approximation it is possible to separate the degrees of 
freedom associated with the relative position r and relative velocity v of the particles from the center of mass degrees 
of freedom. For a system of particles with different masses this approximation is only possible up to the c~ 2 order, 
whereas for a system with equal charge- to-mass ratio (with Z\e/mi = Ziejmi ) the approximation goes up to order 
c~ 4 . Using Lagrange's equation of motion, it is then straightforward to obtain a numerical result for the deflection 
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FIGURE 1. (a) Upper panel: Real part (solid line) ofthe orbital integral I(EIA) for y= 1.1 (£j ab ~ 100 MeV/nucleon), calculated 
with a theory containing both retardation and relativistic corrections to the Coulomb trajectory. The approximation used in ref. [5], 
neglecting relativistic corrections to the Coulomb trajectory, is shown by the dashed line. Lower panel: Same plot, but for the orbital 
integral I(E2, 2). (b) Percentage correction in the calculation of the cross sections for Coulomb excitation of the 0.89 MeV state in 
40 S + 197 Au collisions as a function of the bombarding energy. The solid line (NR) corresponds to the use of the non-relativistic 
orbital integrals [4] compared to the "exact" calculation of ref. [9]. The same is plotted for the other two cases: (R) with retardation 
effects only [5], and (RR) with the retardation effects plus an approximate recoil correction [5]. For more details, see ref. [5]. 



angle and the elastic cross section. Explicit expressions for these Lagrangians are given in ref. [2]. The authors show 
that the scattering angle increases by up to 6% when relativistic corrections are included in 208 Pb + 208 Pb collisions 
at 100 MeV/nucleon. The effect on the elastic scattering cross section is even more drastic: ~ 13% for center-of-mass 
scattering angles around 0-4 degrees. 

Another result obtained in ref. [2] is that the effects of relativity can be taken into account in a much simpler way 
when the projectile is a light particle scattering on a heavy target (e.g., 1 'Li + 208 Pb). In this case, the main contribution 
of relativistic corrections (up to 80% of the total effect) is due to the changes in masses, an easily implementable 
correction. Only about 20% of the relativistic effects are due to magnetic interactions and retardation. Then, a good 
approximation for the elastic scattering is given by 
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with 77 = f5Lc/ZiZ2e 2 , L being the angular momentum of the system, and j3 = v/c. The function is always 
positive for the relevant scattering angles. Therefore one concludes that the relativistic corrections in the elastic 
Coulomb cross section are always negative if they are parametrized as a function of the initial velocity. On the other 
hand if one fixes the initial kinetic energy, instead of the velocity, the corrections will be always positive, as shown in 
ref. [2]. 

The study mentioned above is useful for an analysis of experimental data on inelastic nucleus-nucleus scattering in 
intermediate energy collisions. The extension of these calculations to all orders in j3 is an intriguing problem and an 
appropriate tool could be the use of the action-at-a-distance electrodynamics of Fokker, Wheeler and Feynman [3]. 
Certainly, more studies in this direction are needed and the present focus on nuclear reactions at intermediate energies 
seem to be a good opportunity to investigate such effects. 



3. ANOTHER TRACTABLE CASE: COULOMB EXCITATION 



The theory of Coulomb excitation in low-energy collisions is very well understood [4]. It has been used and improved 
for over thirty years to infer electromagnetic properties of nuclei and has also been tested in experiments to a high 
level of precision. A large number of small corrections are now well known in the theory and are necessary in order to 
analyze experiments on multiple excitation and reorientation effects. 

In the case of relativistic heavy ion collisions pure Coulomb excitation may be distinguished from the nuclear 
reactions by demanding extreme forward scattering or avoiding the collisions in which violent reactions take place. The 
Coulomb excitation of relativistic heavy ions is thus characterized by straight-line trajectories with impact parameter 
b larger than the sum of the radii of the two colliding nuclei. A detailed calculation of relativistic electromagnetic 
excitation on this basis was performed by Winther and Alder [5]. As in the non-relativistic case, they showed how one 
can separate the contributions of the several electric (EX) and magnetic (MX) multipolarities to the excitation. Later, it 
was shown that a quantum theory for relativistic Coulomb excitation leads to minor modifications of the semiclassical 
results [6]. In Refs. [7, 8] the inclusion of relativistic effects in semiclassical and quantum formulations of Coulomb 
excitation was fully clarified. 

Recently, the importance of relativistic effects in Coulomb excitation of a projectile by a target with charge Z2, 
followed by gamma-decay in nuclear reactions at intermediate energies was studied in details. The Coulomb excitation 
cross section is given by 
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where B(nX,Ii — > //) is the reduced transition probability of the projectile nucleus, nX = El, E2, Ml,... is the 
multipolarity of the excitation, and fi = —A, — A + l,...,X. The orbital integrals S(7tX,jX) contain the information 
about relativistic corrections on the relative motion between the nuclei, as well as the relativistic effects on the 
excitation mechanism (e.g. retarded Coulomb interaction). Inclusion of absorption effects in S(nX,}i) due to the 
imaginary part of an optical nucleus-nucleus potential where worked out in ref. [8]. These orbital integrals depend on 
the Lorentz factor 7 = (1 — v 2 /c 2 )~ l l 2 , with c being the speed of light, on the multipolarity %X\l, and on the adiabacity 
parameter E, (b) = (Ofib/ yv < 1, where (Ofi = (Ef — £,) /h is the excitation energy (in units of h) and b is the impact 

parameter. It is often more convenient to write the orbital integrals in another form, e.g. S(EX,jJ.) = —^j- I{EX,jl), 

va A 

where is a geometrical factor, depending only on X and jii, a = Z^e 2 / 'ymv 2 \ and I(EX,lx) is now the orbital 
integral in an dimensionless form. 

Figure 1(a) (from [9]) shows the effect of relativity on the Coulomb interaction and classical trajectory corrections 
in nuclear collisions at intermediate collisions. The comparison is made in terms of the variable £ which is the 
appropriate variable for high energy collisions. Only for ^ < 1 the expressions in the relativistic limit reproduce the 
correct behavior of the orbital integrals. 

TABLE 1. Coulomb excitation cross sections of the first excited state in 38 > 40 > 42 S and 44 > 46 Ar projectiles at 10, 50 
100 and 500 MeV/nucleon incident on gold targets. The numbers inside parenthesis and brackets were obtained with 
pure non-relativistic and straight-line relativistic calculations, respectively. The numbers at the center are obtained 
with the "full" account of relativistic effects, as explained in ref. [9]. 



Nucleus 


[MeV] 


B(E2) 
[e 2 fm 4 ] 


10 MeV/A 50 MeV/A 
Ofc [mb] a c [mb] 


100 MeV/A 500 MeV/A 
Ofc [mb] o c [mb] 


38 S 


1.29 


235 


(492) 500 [651] | (80.9) 91.7 [117] 


(40.5)50.1 [57.1] | (9.8) 16.2 [16.3] 


40 S 


0.89 


334 


(877) 883 [1015] | (145.3) 162 [183] 


(76.1) 85.5 [93.4] | (9.5)20.9 [21.] 


42 S 


0.89 


397 


(903) 908 [1235] | (142.7) 158 [175] 


(65.1) 80.1 [89.4] | (9.9) 23.2 [23.4] 


44 Ar 


1.14 


345 


(747) 752 [985] | (133) 141 [164] 


(63.3)71.7 [80.5] | (8.6) 17.5 [17.6] 


46 Ar 


1.55 


196 


(404) 408 [521] | (65.8) 74.4 [88.5] 


(30.2) 37.4 [41.7] | (5.72) 10.8 [11] 



Table 1 shows the effects of relativistic corrections in the collision of the radioactive nuclei 38,40 ■ S and 44 > 46 Ar on 
gold targets. These reactions have been studied at E[ a b ~ 40 MeV/nucleon at the MSU facility [10]. Table 1 shows the 
Coulomb excitation cross sections of the first excited state in each nucleus as a function of the bombarding energy per 
nucleon. The cross sections are given in milibarns. The numbers inside parenthesis and brackets were obtained with 



pure non-relativistic [4] and relativistic calculations [5], respectively. The minimum impact parameter is chosen so 
that the distance of closest approach corresponds to the sum of the nuclear radii in a collision following a Rutherford 
trajectory. One observes that at 10 MeV/nucleon the relativistic corrections are important only at the level of 1%. At 
500 MeV/nucleon, the correct treatment of the recoil corrections is relevant on the level of 1 %. Thus the non-relativistic 
treatment of Coulomb excitation [4] can be safely used for energies below about 10 MeV/nucleon and the relativistic 
treatment with a straight-line trajectory [5] is adequate above about 500 MeV/nucleon. However at energies around 50 
to 100 MeV/nucleon, accelerator energies common to most radioactive beam facilities (MSU, RIKEN, GSI, GANIL), 
it is very important to use a correct treatment of recoil and relativistic effects, both kinematically and dynamically. At 
these energies, the corrections can add up to 50%. These effects were also shown in Ref. [7] for the case of excitation 
of giant resonances in collisions at intermediate energies. 

We conclude that a reliable extraction useful nuclear properties, like the electromagnetic response (B(E2)-values, 
y-ray angular distribution, etc.) from Coulomb excitation experiments at intermediate energies requires a proper 
treatment of special relativity. The effect is highly non-linear, i.e. a 10% increase in the velocity might lead to a 
50% increase (or decrease) of certain physical quantities. 



4. STRONG INTERACTION: GLAUBER, DWBA AND SEMICLASSICAL METHODS 

The treatment of the strong interaction in nucleus-nucleus collisions at intermediate and high energies is evidently 
much more complicated than the case of Coulomb excitation, described in the previous section. Fortunately, many 
direct nuclear processes, e.g. nucleon knockout, or stripping, elastic breakup (diffraction dissociation), etc, are possible 
to study using the optical limit of the Glauber theory, in which the nuclear ground-state densities and the nucleon- 
nucleon total cross sections are the main input. In fact, this method has become one of the main tools in the study 
of nuclei far from stability [11]. The reason is that the eikonal (or Glauber) methods only use the dependence of 
the scattering matrices, S(b), on the transverse direction, b. Transverse directions are always Lorentz invariants. The 
reason for using S(b), instead of 5(r), traces back to the eikonal scattering wavefunction at the asymptotic region 
(r — oo), 

y Vscatt=S(b) exp(ik-r) , (4) 

where k is the particle's momentum, and r its position. Obviously, the plane wave part of eq. 4 is Lorentz invariant. The 
S-matrix in eq. 4 is S = e\p{ix(b)}, where %(b) is the eikonal phase-shift, given in terms of the interaction potential 
Vby 

x =-hL dzV ^- (5) 

Under a Lorentz transformation to the target system, the coordinate z transforms as z — > JZ- Thus, strictly speaking, the 
S-matrix S(b) is Lorentz invariant only if V transforms as the time-component of a four- vector, i.e. V (r) — ► yV (b, yz) ■ 
The relativistic property described above is most easily seen within a folding potential model for a nucleon-nucleus 
collision: 

V(r) = Jdr' 3 p T (r') v NN (r-r'), (6) 

where pj (r') is the nuclear density of the target. In the frame of reference of the projectile, the density of the target 
looks contracted and particle number conservation leads to the relativistic modification of eq. 6 so that p T (r') — > 
YPt { r 'j_> Yz'), where is the transverse component of r'. But the number of nucleons as seen by the target (or 
projectile) per unit area remains the same. In other words, a change of variables z" = yz 1 in the integral of eq. 6 
seems to restore the same eq. 6. However, this change of variables also modifies the nucleon-nucleon interaction vnn- 
Thus, relativity introduces non-trivial effects in a potential model description of nucleus-nucleus scattering at high 
energies. Colloquially speaking, nucleus-nucleus scattering at high energies is not simply an incoherent sequence of 
nucleon-nucleon collisions. Since the nucleons are confined within a box (inside the nucleus), Lorentz contraction 
induces a collective effect: in the extreme limit y — > °o all nucleons would interact at once with the projectile. This is 
often neglected in pure geometrical (Glauber model) description of nucleus-nucleus collisions at high energies, as it is 
assumed that the nucleons inside "firetubes" scatter independently. 

Assuming that the nucleon-nucleon interaction is of very short range so that the approximation vnn (r — r') = 
Jo 8 (r — r') can be used, one sees from eq. 6 that V (r), the interaction that a nucleon in the projectile has with the 
target nucleus, also has similar transformation properties as the density: V (r) — ► yV (r±,yz), i.e. V (r) transforms as the 
time-component of a four- vector. In this situation, the Lorentz contraction has no effect whatsoever in the diffraction 



dissociation amplitudes, described in the previous sections within the eikonal approximation. This is because a change 
of variables z' = yz in the eikonal phases leads to the same result as in the non-relativistic case, as can be easily checked 
from eq. 5. Of course, the delta-function approximation for the nucleon-nucleon interaction means that nucleons will 
scatter at once, and Lorentz contraction does not introduce any additional collective effect. This is not the case for 
realistic interactions with finite range and collisions at intermediate energies. 

Using the eikonal approach to account for scattering of particles within the projectile, one obtains the wavefunctions 
for the initial and final states as 

% = <h (r ) exp (ik • R) , <F, = 0/ (r) S (b) exp (ik • R) , (7) 

where 0;,/(r) are the initial and final probability amplitudes (wavefunctions) that a particle in the projectile is at a 
distance r from its center of mass. The particle's 5-matrix, S (b), accounts for the distortion due to the interaction. If we 
now assume that the projectile is a two-body system (e.g. a core+valence particle), we get the diffraction dissociation 
formula as follows. The wavefunction of a two-body projectile in the initial and final states is given by 

% = <j>i (r) exp [i (k c • r c + k v • r v )] , V f = <j> f (r) S c {b c ) S v (b v ) exp [i (k' c ■ r c + k' v -r v )] , (8) 

where now j (r) are the initial and final intrinsic wavefunctions of the (core+valence particle) as a function of 
r = Y\ — Y2- The relation between the intrinsic, r, and center of mass, R, coordinates is given in terms of the mass ratios 
j3, = nti/rnp. Explicitly, r v = R + j3 c r and r c = R — j3 v r. The core and valence particle 5-matrices, S c (b c ) and S v (b v ), 
account for the distortion due to the interaction with the target. 

The probability amplitude for diffraction dissociation is the overlap between the two wavefunctions above, i.e. 

^(diff) = j d 3 r c d 3 r v <j)* f (r) (r) 8 (z c + z v ) S c (b c ) S v (b v ) exp [i (q c • r c + q v • r v )] , (9) 

where q c = k'. - k c is the momentum transfer to the core particle, and accordingly for the valence particle. The above 
formula yields the probability amplitude that the projectile starts the collision as a bound state and ends up as two 
separated pieces, in this case, the core and the valence particle (e.g. a proton or a neutron). All the information for 
the dissociation mechanism comes from the knowledge of the 5-matrices, S c and 5 V . The delta-function 8 (Z) was 
introduced in the eq. 9 to account for the fact that the 5-matrices calculated in the eikonal approximation only depend 
on the transverse direction. 

In the weak interaction limit, or perturbative limit, the phase-shifts are very small so that 

S c {b c ) S v (b v ) = exp [i (Xc + % v )\ ^ 1 + iXc + iXv 

= 1 ~~k.f VcT (rt) dZc 'hi VvT (rv) dZv - (10) 

The factor 1 does not contribute to the breakup. Thus, inserting the result above in eq. 9, we obtain 

Vwba) ~^h~ v \ d3r c d3r v 0/ (r) 0/ (r) [V cT (r c ) + V vT (r v )] exp [i (q c • r c + q v • r v )] , (11) 

where the integrals over z c and z v in eq. 10 were absorbed back to the integrals over r t and r v after use of the delta- 
function 8 (z c + z v )- The above equation is nothing more than the plane-wave Born-approximation (PWBA) amplitude. 
However, absorption is not treated properly. For small values of r c and r v the phase-shifts are not small and the 
approximation used in eq. 10 fails. A better approximation is to assume that for small distances, where absorption is 
important, S c (b c )S v (b v ) ~ S(b), where the right-hand side is the 5-matrix for the projectile scattering as a whole on 
the target. Using the coordinates r and R , and defining f/ (ni (r,R) = V c t (r c ) + V„t (r„), one gets for the T-matrix 

?(DWBA) = »»v-A(dwba) J # rdr 'R $* f (r)exp[iq-r}fr(r)U in Ar,R)S(b)exp[iQ-R} . (12) 

In elastic scattering, or excitation of collective modes (e.g. giant resonances), the momentum transfer to the intrinsic 
coordinates can be neglected and the equation above can be written as 



Vvba) = (* H (R) 0c (r) |£Mr,R)| X {+) (R) <h «) , 



(13) 



which has the known form of the DWBA T-matrix. The scattering phase space now only depends on the center of mass 
momentum transfer Q. When the center of mass scattering waves are represented by eikonal wavefunctions, one has 



z (-)*(R) z (+)(R)~S(i)exp[iQ.R] . (14) 

This shows that the PWBA and the DWBA are perturbative expansions of the diffraction dissociation formula 9. 

In DWBA (or in the eikonal approximation, eq. 14), b does not have the classical meaning of an impact parameter. 
To obtain the semiclassical limit one goes one step further. By using eq. 12 and assuming that R depends on time so 

that R = (b,Z = vt), the semiclassical scattering amplitude is given by A^?^ (b) = f d 2 b a[^£ lass) 0) exp(/Q • b), 
where 

„('-/) 



°(semiclass) 



(b) = i S(b) J dtd\ exp(^ / t) f f (r)U int (r,t)^(r) , (15) 

where Q Z Z = (Off t was used. 

The semiclassical probability for the transition (i — > /) is obtained from the above equations after integration over 



2 

^ (semiclass) 

(b) 



, with b having now the explicit meaning of an impact parameter. Thus, 

a (lemiciass) (^) * s tne semiclassical excitation amplitude. Equation 15 is well-known (for example in Coulomb excitation 
at low energies, where Ui„ t = Uc) except that the factor S (b) is usually set to one. In high energy collisions it is crucial 
to keep this factor, as it accounts for refraction and absorption at small impact parameters: \S (b) \ 2 = exp 2;^ lmag ' , 

where j( ima s) is calculated with the imaginary part of the optical potential. The derivation of the DWBA and 
semiclassical limits of eikonal methods can be easily extended to higher-orders in the perturbation V. The eikonal 
method includes all terms of the perturbation series in the sudden-collision limit. 

The developments presented in this section show that the DWBA calculations of nuclear excitation, and the higher 
order terms, are implicitly included in the eikonal models. However, there is a subtle link to the optical potential, which 
makes the theory Lorentz covariant. Usually, the optical limit of the Glauber-eikonal series is used. In this model, no 
explicit reference to a nuclear potential is done; only the nucleon-nucleon cross sections and nuclear densities are 
used as input (see, e.g. ref. [12]). As shown above, this is also not a guarantee of Lorentz covariance. Even worse is 
the fact that often DWBA (and higher-order, e.g. CDCC [13, 14, 15]) calculations are used without consideration of 
relativistic effects. In the next section I give an example of a method (relativistic CDCC) which incorporates relativistic 
corrections in a continuum-discretized basis [16]. 



5. RELATIVISTIC CONTINUUM DISCRETIZED COUPLED-CHANNELS 

Let us consider the Klein-Gordon (KG) equation with a potential Vo which transforms as the time-like component of 
a four-vector [17] (here I use the notation h = c=\). For a system with total energy E (including the rest mass M), 
the KG equation can be cast into the form of a Schrodinger equation (with h = c = 1), (V 2 + k 2 — t/) *P = 0, where 
k 2 = (E 2 -M 2 ) and U = V {2E - V ). When V <M, and E ~M, one gets U = 2MV , as in the non-relativistic case. 
The condition Vo <C M is met in peripheral collisions between nuclei at all collision energies. Thus, one can always 
write U = 2EVq. A further simplification is to assume that the center of mass motion of the incoming projectile and 
outgoing fragments is only weakly modulated by the potential Vo. To get the dynamical equations, one discretizes the 
wavefunction in terms of the longitudinal center-of-mass momentum k z , using the ansatz 

* = (z,b) exp(ft oZ ) h a (<?) • (16) 

a 

In this equation, (z,b) is the projectile's center-of-mass coordinate, with b equal to the transverse coordinate. <j> (|) 
is the projectile intrinsic wavefunction and (k,K) is the projectile's center-of mass momentum with longitudinal 
momentum k and transverse momentum K. There are hidden, uncomfortable, assumptions in eq. 16. The separation 
between the center of mass and intrinsic coordinates is not permissible under strict relativistic treatments. For high 
energy collisions we can at best justify eq. 16 for the scattering of light projectiles on heavy targets. Eq. 16 is only 
reasonable if the projectile and target closely maintain their integrity during the collision, as in the case of very 
peripheral collisions. 

Neglecting the internal structure means <j) ka (§) = 1 and the sum in eq. 16 reduces to a single term with a = 0, 
the projectile remaining in its ground-state. It is straightforward to show that inserting eq. 16 in the KG equation 



(V 2 + k 2 - 2EV ) ¥ = 0, and neglecting V 2 y Q (z, b) relative to ikdz^b (z, b), one gets i/tofe^o (z, b) = £Vo^o {z, b), 
which leads to the center of mass scattering solution o5^o fob) = exp [— iv~ 

1 SL,dn! V Q (z',b)] , with v = k/E. Using 
this result in the Lippmann-Sch winger equation, one gets the familiar result for the eikonal elastic scattering amplitude, 
i.e. /o = —i(k/2n) / db exp(z'Q-b) {exp[;'j(b)] - 1}, where the eikonal phase is given by exp[/j(b)] = ^b(°°,b), 
and Q = K 7 K is the transverse momentum transfer. Therefore, the elastic scattering amplitude in the eikonal 
approximation has the same form as that derived from the Schrodinger equation in the non-relativistic case. 

For inelastic collisions we insert eq. 16 in the KG equation and use the orthogonality of the intrinsic wavefunctions 
(j>k a (§). This leads to a set of coupled-channels equations for y a - 

(v 2 +k 2 ) y a e kaZ = £(« \u\ a 1 ) y a , e k »'\ (17) 

a 

with the notation \a) = \tyk a )- Neglecting terms of the form V 2 y a (z,b) relative to ikdz^a (z,b), eq. 17 reduces to 



iv d.y a (z,b) ^(a\Vo\a>) y a ,{zM A'*-**)'. (18) 

" z a 1 

The scattering amplitude for the transition — > a is given by 

fa (Q) = ~ j db exp (jQ • b) [S a (b) - 8 afi ] , (19) 

with S a (b) = y a (z = °°,b). The set of equations 18 and 19 are the relativistic-CDCC equations (RCDCC). 

The RCDCC equations have been used [16] to study the dissociation of 8 B projectiles at high energies. The energies 
transferred to the projectile are small, so that the wavefunctions can be treated non-relativistically in the projectile 
frame of reference. In this frame the wavefunctions are described in spherical coordinates, i.e. \a) — \jlJM), where 
j, I, J and M denote the angular momentum numbers characterizing the projectile state. Eq. 18 is Lorentz invariant if 
the potential Vo transforms as the time-like component of a four-vector. The matrix element (a \Vq\ a') is also Lorentz 
invariant, and one can therefore calculate them in the projectile frame. 

The longitudinal wavenumber k a ~ (E 2 —M 2 ) 1 ! 2 also defines how much energy is gone into projectile excitation, 
since for small energy and momentum transfers k' a —k a = (E' a — E a )/v. In this limit, eqs. 18 and 19 reduce to 
semiclassical coupled-channels equations, if one uses z = vt for a projectile moving along a straight-line classical 
trajectory, and changing to the notation y a (z,b) — a a (t,b), where a a (t,b) is the time-dependent excitation amplitude 
for a collision wit impact parameter b (see eqs. 41 and 76 of ref. [18]). The full version of eq. 19 was used in ref. [16], 
with relativistic corrections in both the Coulomb and nuclear potentials. 

If the state \a) is in the continuum (positive proton+ 7 Be energy) the wavefunction is discretized according to 
\tt\E a ) = S dE' a T(E a ) \a;E' a ), where the functions T(E a ) are assumed to be strongly peaked around the energy 
E a with width AE. For convenience the histogram set (eq. 3.6 of ref. [19]) is chosen. The inelastic cross section is 
obtained by solving the RCDCC equations and using da/dQ.dE a = \ f a (Q)| 2 T 2 {E a ). 

Figure 2 shows the relative energy spectrum between the proton and the 7 Be after the breakup of 8 B on lead targets 
at 83 MeV/nucleon. The data are from ref. [20]. In this case, the calculation was restricted to b > 30 fm. The dotted 
curve is the first-order perturbation calculation, the solid curve is the RCDCC calculation, and the dashed curve is 
obtained with the replacement of 7 by unity in the nuclear and Coulomb potentials. The difference between the solid 
and the dashed-curve is of the order of 4-9%. 



6. CONCLUSIONS 

The consequence of neglecting relativity in nuclear reactions at intermediate energies is not easy to access. The 
inclusion of relativity introduces non-trivial effects in semiclassical, DWBA, eikonal, and continuum discretized 
coupled-channels calculations. Nuclear collisions are up to now the most used probe of the internal structure of rare 
nuclear isotopes. To my knowledge, most experiments have been analyzed using non-relativistic theoretical methods. 
It might be necessary to review the results of some of these data, using a proper treatment of the relativistic corrections 
in the theoretical calculations used in the experimental analysis. Other improvements of the formalisms presented here 




and elsewhere needs to be assessed. The relativistic effects in the nuclear interaction has also to be studied in more 
depth. 

Special relativity [1], one of the most precious theories of Eintein's legacy, still remains a source of intriguing 
effects, not always easy to tackle. Nuclear physics is full of such examples. 
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